Normal and anomalous diffusion in random potential landscapes 
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A relation between the effective diffusion coefficient in a lattice with random site energies and 
random trasition rates and the macroscopic conductivity in a random resistor network allows for 
elucidating possible sources of anomalous diffusion in random potential models. We show that sub- 
diffusion is only possible either if the mean Boltzmann factor in the corresponding potential diverges 
or if the percolation concentration in the system is equal to unity (or both), and that superdiffusion 
is impossible in our system under any condition. We show also other useful applications of this 
relation. 

PACS numbers: 05.45.-a; 05.60.Cd 
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A classical particle's diffusion in a random potential or 
hopping on a lattice with disordered site energies E^ is 
a versatile theoretical model with a wide range of appli- 
cations [lj. The particle's motion typically corresponds 
to normal diffusion, but can get subdiffusive in the pres- 
ence of deep traps or in the case of infinte contrast when 
approaching percolation transition. A question arises, 
whether there can be other cases leading to sub-diffusion 
except for these two (or combinations thereof). In what 
follows we show that the diffusion coefficient in a discrete 
disordered lattice is always finite (i.e. that no superdiffu- 
sion can be observed) but may vanish (possibly giving rise 
to subdiffusion) . Independently of the particular distri- 
bution of nonzero transition rates, this is possible either 
if the percolation threshold in the corresponding network 
is unity (e.g. in one dimension, on finitely ramified frac- 
tals, or when we are already at percolation threshold) or if 
the mean Boltzmann factor (exp(— Ei / kT)) diverges (or 
both). According to the Arrhenius law, the last situation 
corresponds to the divergence of the mean sojourn time 
at a site, pertinent to trapping. 

We start from the Master equation for the probabilities 
Pi to find a particle at a site i of a lattice (with lattice 
spacing a) 
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(■WijPj - WjiPi) 



(1) 



where Wij are transition rates from site j to site i differ- 
ent from zero only for nearest neighbors. Eq. (|T|) can ei- 
ther follow from some microscopic scheme or be obtained 
by a discretization of the Fokker-Planck equation for the 
overdamped motion in a continuous potential. The sys- 
tem is taken to be homogeneous and isotropic in statisti- 
cal sense. This requirement excludes underdamped cases 
for which the velocities and coordinates enter differently, 
thus leading to anisotropy of the state (phase) space. In 
what follows we consider a d-dimensional lattice with to- 
tal of M ^> 1 sites assigned energies Ei being identically 
distributed random variables. We assume that the sys- 
tem is isothermic and possesses true thermodynamical 
equilibrium under appropriate boundary conditions, i.e. 
that the transition rates fulfill the detailed balance condi- 



tion w^p® = Wjip\ at equilibrium (the superscript will 
denote the corresponding value at equilibrium through- 
out the work). The transition rates < Wij < oo are 
not necessarily bounded from above, and some of them 
may be put to zero to mimic percolation situations. The 
values of p\ are given by p\ oc hi = exp(—Ei/kT), where 
bi denotes the Boltzmann factor, T is the temperature 
and k is the Boltzmann constant. 

Our discussion maps the initial problem onto the one 
for random resistor-capacitor networks. Let be the 
corresponding conductivities of the bonds, and (gij) EM 
be the effective conductivity of a corresponding network 
in the static regime. Then the effective diffusion coeffi- 
cient in a network follows as 
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The statements done in the first paragraph are then 
demonstrated by using the results from the theory of elec- 
tric circuits and from the percolation theory. Some other 
useful applications of Eq. @ are shown. 

Eq.Q by itself is not new, but we give here its physical 
derivation which stresses its general applicability and its 
connection with thermodynamics. Thus, the discussion 
for the case of a barrier model (all Ei are the same, but 
the transition rates fluctuate) is contained in Ref. [lj, 
Eq.(2.15). Moreover, Eq.© naturally appears when ap- 
plying effective medium approximation (EMA), like the 
one of 0. A derivation for a continuous case (Langevin 
description in Ito interpretation) is given in Q . Note that 
the Ito prescription may correspond to the trap model in 
the discrete case i.e. to a situation different from 
the one of Ref. [lj . In our work we confine ourselves to 
a discrete setup which allows for the application of the 
theory of electric circuits for the analysis of the results 
(although we make a limiting transition to continuum to 
illustrate some outcomes of the approach). 

We start by rewriting Eq.([T]) as an equation for mean 
numbers ("concentrations") of non- interacting particles 
at the corresponding sites, hi — (wijUj — Wjifii), con- 
nected with probabilities via rii = Npi with N being 
the total number of particles. In equilibrium all n® are 
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proportional to the Boltzmann factors, = Cbi with 
prefactor C depending on the number of particles, on 
the system's size and on distribution of 6j. Putting the 
detailed balance condition into the form w^n® — w^n® 
we denote Wijn® = <7,j where Qij = gji is now a prop- 
erty of the bond. Using this notation we rewrite Eq.([T]) 
as an equation for the temporal evolution of activities 
Ci — n i/ n i ( see the Appendix): 

Ci = -o tiijtj ~ 9jiCi) ■ (3) 

Eq. ([3]) is formally equivalent to the evolution equation of 
node potentials in a random resistor-capacitor model [l| , 
with conductivities gij and capacitances n®. 

Let us now calculate the effective diffusion coefficient 
provided it exists (i.e. the system homogenizes at large 
scales). For random resistor networks the homogeniza- 
tion of conductivity is mathematically proved for local 
conductivities bounded from above and from below, see 
Q and references therein. The boundness from below 
excludes the conductor-isolator percolation model, but 
homogenization still holds provided the system perco- 
lates [7[. Physically, it is known that the conductor- 
superconductor system homogenizes below the percola- 
tion threshold for superconductor Q. 

We mimic a stationary experiment on measuring the 
diffusion coefficient via the first Fick's law: The system, 
in a form of a membrane of thickness L and of transversal 
dimension W separates two reservoirs, the left one with a 
well-stirred solution of particles at concentration v\ , the 
other one with a slightly lower concentration v r . The 
membrane is considered as impermeable for the solvent, 
and both concentrations are kept constant. The constant 
particles' current / through the membrane is measured 
and connected with the mean diffusion coefficient inside 
it. Since in general a jump of the (free) energy per par- 
ticle can form on a contact between the membrane and 
the solution (e.g. when the fluid is a good solvent for dif- 
fusing particles and the membrane is, on the average, a 
bad one, or other way around) the effective diffusion co- 
efficient inside the membrane has to be defined through 



W^dm) - (n r )) w 

where (ni) and (n r ) are the mean particle concentrations 
in the layers of the membrane in immediate contact with 
the solution, see Fig.l. In a stationary state Ci = 0. 
Moreover, in the thermodynamical limit L — > oo the per- 
meability of the membrane and thus the current tend to 
zero. We will call this situation "quasi-equilibrium" in 
what follows. 

The contact with solution is modelled by additional 
arrays of sites to the left and to the right from the mem- 
brane, with constant particles' concentrations and con- 
stant energies Eq which can be chosen arbitrarily (Eq 




FIG. 1: A schematic illustration of the situation considered 
in the text: the disordered medium in contact with two reser- 
voirs, the mean concentration at different positions and the 
lattice model applied. 



defines the quality of the solvent). These additional 
sites are connected to the ones on the membrane's sides 
via extremely high transition rates fulfilling the detailed- 
balance condition. In this case a local equilibrium be- 
tween the surface sites and the solutions persists inde- 
pendently on the particles' distribution inside the bulk. 
Due to this the activities of the surface sites are all 
equal to Cl = Ai>i and Cr = Av r at the left resp. right 
boundary of the membrane, where the prefactor A de- 
pends on Eq. Thus, in the leftmost layer are pro- 
portional to Qn® — Avin { \ and in the rightmost layer to 
Crrf = Ais r n® so that the mean concentrations in the 
layers are rt; = Avi(rq) and n r = Ah> r (n®}, where we 
assume that the distribution of the site energies in the 
surface layers is the same as in the bulk. We then calcu- 
late the corresponding total current I through the system 
noting that the equations for the currents and activities 
in a stationary state are the same as the ones given by 
the Kirchhoff's laws for an electric circuit. Making such 
a reinterpretation we see that / = g*{W d ^ 1 / L)(Q — Cr) 
where g* = (g) em is the effective conductance (conduc- 
tivity of a bond in the effective ordered medium with the 
same total conductivity as our heterogeneous one) , where 
the subscript EM denotes the effective medium mean. 
Therefore D* = a 2 g*/(n®). The prefactor a 2 is intro- 
duced to restore the dimension as follows from Eq. (|4|), 
when passing from distances L measured in lattice units 
to distances measured in centimeters. We note that since 
n° are proportional to the Boltzmann factors, and since 
rescaling of all g^ by a constant factor leads to chang- 
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ing g* by the same factor, the proportionality factor C 
cancels out; this gives Eq.©. 

The result is rather transparent. If we are able to 
measure the effective conductivity of the system, we 
can connect it with the effective mobility /i* (and thus 
with the diffusion coefficient) via Nernts-Einstein equa- 
tion a* — noqfi* = n^qD* /kT, where hq is the equilib- 
rium concentration of particles with charge q. Reverting 
this expression we get D* oc a* /no, which is essentially 
Eq.®. 

One may argue that the correct way is to define D* 
through the gradient of the coarse-grained concentration, 
and not via the total concentration difference. As we 
show in the Appendix, this definition leads to the same 
result since local concentrations and local activities de- 
couple under quasi-equilibrium (but only under this con- 
dition!). 

Let us first discuss some applications of Eq.© other 
than discussed in the preface. Eq.© gives the possibil- 
ity to obtain the universal bounds on the effective dif- 
fusion coefficient based on those for the effective con- 
ductance, i.e. the universal Wiener bounds Q and the 
tighter Hashin-Shtrikman bounds for isotropic systems 
10], as well as to generalize some exact results for two- 
dimensional systems based on duality . In this presen- 
tation we concentrate on continuum models where Eq.(fTj) 
arises from discretization of a Fokker-Planck equation for 
p(x,t): p = DAp+(D/kT)V(VU(x)p) with constant dif- 
fusion coefficient D and disordered potential U (x) . The 
details of calculations are given in the Appendix. 

The universal Wiener bounds for the conductance are 
given by (g^ 1 )" 1 < (9ij) EM < (9ij), In our cases this 
corresponds to 



a 2 WQ 



(exp 



(Ei/kT))(exp(-Ei/kT)) 



< D* < a 2 w . (5) 



Note that the lower bound reproduces the exact result for 
the one-dimensional system with random potential and 
constant diffusion coefficient. 

In Figs. [5] and [3] we plot the Hashin-Shtrikman 
bounds for two cases: the case of the binary disorder 
Ei = Ei with probability p and Ei = Ei with prob- 
ability 1 — p, and the case of Ei possessing an expo- 
nential distribution with cutoffs, the one with density 
P{Ei) = Pe'^le-P^ - e"* 3 ^]- 1 for E 2 < E t < E x 
and vanishing elsewhere (apart from cutoffs this distri- 
bution is reminiscent of the exponential energy distribu- 
tions leading to CTRWs). As a comparison, the results 
of the effective medium approximation (EMA), see Refs. 
121 Il3j , are shown. The results are plotted as the func- 
tion of a contrast x = exp[(i?2 — E\)/kT\ being the ratio 
of the maximal and the minimal value of gtj . 

Note that the result for the effective diffusion coeffi- 
cient for symmetric binary case in 2d is essentially exact 
since in this case (g) em — \j9a9b due to the duality re- 
lation [13. 
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FIG. 2: EMA result and Hashin-Shtrikman bounds for the 
symmetric binary case (p = 1/2) vs. contrast x for: d — 2 (a) 
and for d — 3 (b). 
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FIG. 3: EMA result and Hashin-Shtrikman bounds for trun- 
cated exponential P(Ei) vs x for d — 2 (a) and d — 3 (b). 



In the limit of very strong disorder D* may vanish or 
diverge. In the first case D = 0, the system either does 
not show any transport (does not percolate) or shows 
anomalous transport slower than diffusion (i.e. shows 
subdiffusiom) . In the second case it might show superdif- 
fusion. 

If D* vanishes, it can do so either because the nu- 
merator (wjiexp(—Ei/kT)) EM vanishes or because the 
denominator (exp(—Ei/kT)) diverges, as well as in the 
cases when both possibilities are realized simultaneously 
(which may give rise to subdiffusion of mixed origins [H). 
We discuss the conditions under which the correspond- 
ing behavior may take place. If D* diverges, it can do so 
because the numerator diverges, or because the denom- 
inator vanishes, or both. As we proceed to show, none 
of these possibilities can be realized. This excludes not 
only superdiffusion, but also the "compensated" cases of 
normal diffusion when both numerator and denominator 
vanish or diverge simultaneously. 

For further discussion we first recapitulate the follow- 
ing properties of percolation systems: (i) The mixture of 
resistors with given finite conductivity (at concentration 
p) with insulating bonds (of zero conductivity) at concen- 
tration 1 — p possesses zero conductance below the perco- 
lation threshold p c and finite conductance above it. The 
corresponding system homogenizes at scales above the 
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correlation length This homogenization also takes 
place for arbitrary distribution of the conductivities of 
the resistors Q- Similarly, (ii) The mixture of resistors 
with given finite conductivity (at concentration p) with 
superconducting bonds (at concentration 1 — p) possesses 
finite conductance below the percolation threshold p s c for 
superconducting bonds, with 1 — p s c = p c , and infinite 
conductance above it. These properties do hold not only 
for the Bernoulli percolation model but also in the case 
when the short-range correlations in the occupation prob- 
abilities of the bond by the corresponding resistors / in- 
sulators / superconductors are present. This statement 
is a (silently assumed) basis of all renormalization group 
approaches in percolation. 

U sing the results of the theory of electric circuits, see 



e.g. [14|, we show that the total conductivity of a resistor 
network is a non-decaying function of the conductivity 
of each particular bond. Let us consider the system as 
placed between two "superconducting" bars considered 
as a terminal 1 of the system. Let us consider the poles 
i and j between which g^ is switched as terminal 2. Us- 
ing the theory of two-terminal circuits we calculate the 
input impedance (total conductivity) z- m as a function 
gij-. z- m = z n - zi22 2 i/(;z22 + gij), where z a p are the 
elements of the impedance matrix of the system. For a 
system of reciprocal passive elements (no batteries, no 
diodes) this matrix is non-negatively definite and sym- 
metric as a consequence of non-negative heat production 
and of reciprocity theorem. In the case of pure resis- 
tor network the matrix is real. Thus, 2:11,2:22 > and 
2122:21 = (212) 2 > 0, so that z- m is a non-decaying func- 
tion of g^ . 

Now we show that the numerator never diverges. We 
fix some q < 1 — p s c = p c and declare the fraction q of 
bonds (starting from the ones with largest g) to be super- 
conductive. The lowest conductivity of a changed bond 
is <7min- The superconducting bonds are non-percolating 
by construction, and the conductance of the remaining 
system is finite, being smaller than a conductance of the 
resistor-superconductor mixture where all conductivities 
are put to g m - ln . Thus the numerator can only diverge if 
p c = 0, i.e. never in finite dimension. 

The numerator does not vanish for a system with per- 
colation concentration p c < 1. Let us remove a portion 
q < 1 — p c of bonds with smallest g without destroying 
percolation and denote the largest removed conductivity 
by jmax > 0. The rest of the system percolates and has 
a conductance which is larger then the conductance of a 
two-phase system constructed of resistors with g — g max 
and g — 0, which is nonzero since we are above perco- 
lation threshold. Thus the numerator can only vanish if 
p c = 1, and no bonds can be removed. 

The denominator in Eq.([2]) can diverge if the corre- 
sponding mean value of the Boltzmann factor diverges. 
Since hi = exp(— Ei/kT) is proportional to the sojourn 
time at a site i in equilibrium, this corresponds to di- 



verging mean sojourn time at a site, i.e. to a trap model, 
which in high dimensions is equivalent to CTRW with a 
broad distribution of waiting times. 

The denominator can not vanish. Let p(E) be the 
probability density of Ei, and Em > —00 its median, 
J E ™p(E)dE = 1/2. Since (6) = (exp(- E t / 'kT)) = 
~ E / kT p(E)dE where the integrand is non- 
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negative, we have (b) > J E ^ e- E / kT p{E)dE > 

exp{-E M /kT) ff£ P (E)dE = (1/2) cxp(-S M /fcT) 
since exp(—E/kT) is monotonically decaying. 

Summarizing our findings we state that there exists an 
exact correspondence between the effective diffusion co- 
efficient in a random potential and macroscopic conduc- 
tivity in a random resistor model. This simple relation 
allows us to obtain exact bounds on the effective diffusion 
coefficient. It also allows for elucidating possible sources 
of anomalous diffusion in such model. Thus, the subd- 
iffusion is possible either if the mean Boltzmann factor 
of the corresponding potential diverges (energetic disor- 
der) or if the percolation concentration in a system is 
equal to unity, i.e. if the system is already at the perco- 
lation threshold, in one dimension, or on finitely ramified 
fractals (structural disorder), and that superdiffusion is 
impossible in our system under any condition. 

Acknowledgement. The work was supported by 
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Appendix 
A - Activities and chemical potentials 

The idea behind the approach is based on the fact that 
the probabilities pi depend on the system's size M and 
therefore are awkward to use in mimicking a macroscopic 
experiment. Thus, according to normalization condition, 
J2itiPi = 1> an d in equilibrium p° t = &i/X^_ x &i or 
pf = Af _1 6j/(6i) provided the corresponding mean ex- 
ists. On the other hand, rii = Npi and n° are the "in- 
tensive" variables, e.g. n® — (N j M)bi / (hi) , and do not 
depend on the system's size provided the mean particle 
concentration N/M is kept constant. 

To see that £i introduced in the main text are indeed 
activities one can proceed as follows. Adding a weak ex- 
ternal potential (e.g. electric field giving rise additional 
potential energy Vj at site j) changes the equilibrium con- 
centrations and hence the corresponding equation to hi = 
Ej L>« eMVj/kT) - gjiim/n?) exp(U 4 /fcT)] . 

Here Q = (rii/n®) exp(Vi/kT) can indeed be interpreted 
as activities, since /i; = kT hi Q = fcTln^- + Vi have a 
form of chemical potentials on the sites. 

The master equation can be interpreted as a combi- 
nation of the local continuity equation hi = Jij and 
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the local linear response equation Jy = gij(Q ~ Ci) with 
Cj = njexp(Vj/kT)/rij. Close to equilibrium and for 
very small potentials V the values of ( are close to unity, 
so that the deviations of chemical potentials from their 
equilibrium values of zero are small. Therefore, close to 
equilibrium, 



and 



J i 



9v_ 
kT 



(pj -(h)- 



The external potential Vj is set to zero in the main text. 

It is interesting to note that not the chemical poten- 
tials but the activities at the sites are reinterpreted as 
electric potentials, which leads to differences in the be- 
havior of the random resistor-capacitor model and the 
random potential model far from equilibrium, i.e. for 
large concentration gradients. 



B - Decoupling of local concentrations and local 
activities under quasiequilibrium 

In a stationary situation corresponding to quasiequi- 
librium we have 



~ Ci = 0, 



(6) 



where the sum runs over the nearest neighbors of the site 
i. Moreover, in the thermodynamical limit L — > oo the 
permeability of the membrane tends to zero, and all local 
currents in it as well: 



Ji 



ffy(O-Ci)-K). 



(7) 



Now we show that in quasi-equilibrium the local con- 
centrations and the local activities ("potentials") tend 
to be independent from each other (although the lo- 
cal transition rates are correlated with the Boltzmann 
factors and thus local concentrations). To see this we 
note that according to Eq.© the activity (potential) 
at a site i is the weighted arithmetic mean of the ac- 
tivities of the neighbors it is connected to, and Eq.((7J) 
states that the difference between the activities of the 
connected sites under quasi-equilibrium tends arbitrarily 
small. Therefore our system can be considered as com- 
posed of large, practically equipotential regions whose 
potential hardly fluctuates around its mean depending 
on the region's position. In these regions the concentra- 
tions tii — dri® can be averaged over the physically small 
volume still containing the large number of sites, so that 
(n) = C*{n^). We note that this last property does not 
rely on homogenization or on the isotropy and will hold 



even if our system is built of independent parallel or in- 
terwoven wires! 

Let us now assume that at large scales the correspond- 
ing electric system homogenizes. In this case the total 
voltage profile (which in our case is mapped on the pro- 
file of activities) obtained by the coarse graining (moving 
average) of the voltages (activities) over macroscopic do- 
mains of the system which are large enough compared 
to the lattice spacing but small compared to L follows 
a linear behavior (similar to those of the effective con- 
centration in Fig.l). Then the coarse-grained activity 
is a linear function of the coordinate, and thus also the 
coarse grained concentration is linear in coordinate, with 
the proportionality factor (n°) between the both. 



C - Discretization of a continuous Fokker-Planck 
equation with random potential 

We start this section by showing the discretization pro- 
cedure linking the Fokker-Planck equation to the master 
equation. The equation in a continuous d-dimensional 
space is discretized using a regular hypercubic lattice 
with lattice constant a considered sufficiently small. The 
sites are identified by the vectors i, the site energies E\ 
are given by corresponding values of the random poten- 
tial U(x) at the position of site i, x = i = (ii, ...,id) — 
a(ni, rid) where n& is an integer number and k runs 
from 1 to d . Correspondingly, the derivatives dk along a 
certain direction k are linked to the forward differences 
A k in the following way: 

d k p(x) -> -A^ = -[p i+ak - Pi] 

a a + 

d 2 kP ( X ) -> -Al Pi = ^[p i+ak ^2p i +p i _ ak ] (8) 

where k is the unit vector in the fe-th direction. Hence, 
from 



p(x) 



DV 2 p(x) + D/3V[VE/(x) -p(x)] = 

d 

[fl£p(x) + (3d k [dkU( X ) ■ p(x)] 



fc=l 



follows 



d 



fc=l 



~^((^+«Jt-^)pl 



(9) 



We consider the potential to be a slowly varying function 
of the position; this allows to substitute the remaining 
partial derivative with half of a double step forward dif- 
ference, the consequent error being of the smaller order 
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of magnitude &? k E. 



D 



Pi = ~2 



a 

E [fo+ak- 2 Pi+Pi-ak} 



(10) 



fc=l 



(( £; i+2afe ^ ^i+afe^i+afe _ (^i _ E i-ak)Pi-ak 



Neglecting again second order terms we replace the quan- 
tity ( E i+2ak - E i+ak) With ( E i+ak ~ S 0. and add and 

subtract the term {E i+a f, — Ei)pi 



P» = ^ S |W«fc - 2 ^ + Pi-ak\ + 



k=l 



+ f ((^i+afe - ^i)Pi+afe - - E i-ak)Pi-ak + 

+ ( E i+ak E i)Pi ( E i+ak E i)p)] = (11) 

^E [bi+afc- 2 Pi+Pi-afe] + 
fe=l 

+ f ((^i+afe - ^i)Pi+afe - - E i- ak )Pi-a k + 

+( E i+ak E i)Pi + - = (12) 

5e[^( i -§(^-^))+ 



fc=i 



+». i (l--{E i -E. i 



(13) 



Pi(2-f(^, +0 ft-^)-f(^-a£-^i)) 



0, 



Since the difference between the energies at neighbouring 
sites is small, in the lowest nonvanishing order we can put 

D d 

* = ^E[[Pi+afe e " l(Si " Bl+ ^ + 
fc=l 

+ Pi- ak e-^ Ei - E >-^ + 

- Pi ( e -K^ + a S -^) +e -|(^-^-^))] (14) 

At this point, the discretized Fokker-Planck equation is 
shown to be equivalent to the master equation 

d 

P 1 = J2[ W i,i+akPi+a k + W i,i-a k Pi-ak + 



fc=l 



- Pi( W i+ak,i + W i-ak,i) 



(15) 



by setting 



w lJ= B-$(E l -E i)=Woe -$(E l -E i ) 

implying that, up to constant factors 



In what follows ^{E\ + can be changed for U(x)/kT. 
Thus, returning to the continuous limit when calculat- 
ing the macroscopic (effective medium) conductance of 
the continuous disordered medium we can take its lo- 
cal conductivity to be g(x) = g cxp(~U(x)/kT) with 
go = D/a 2 . 

D - Calculations pertinent to Figs. 2 and 3 

We now proceed to compute the Hashin-Shtrikman 
bounds for the effective diffusion constant by consider- 
ing two different distributions of the energy variable E\. 



Binary distribution 



In the binary case 



Ei with probability pi 
Ei with probability p2 



P 

1-p 



(16) 



with E 2 < E\. The local conductance follows again a 
binary distribution 



9i,i 



91 oc g e 

92 oc g e 



-PE 2 



with probability Pi = p 
with probability p 2 = 1 — p 

(17) 



or, equivalently, 

p(g) = PiS{g - gi) + P2% - g 2 )- 



(18) 



Here we introduce the following notation for different 
averages which will repeatedly appear in what follows 
(weighted arithmetic, arithmetic and geometric mean): 

(.9) = + 



(9) A = 2(91 +92) 
(9)G = V9192 



(19) 



The notation (f(E)) will be continuously used for the 
mean of any physical quantity / over the energy distri- 
bution. The Effective Medium Approximation then gives 
the effective conductance g* as the solution for the self- 
consistency condition 



(d - 1)5* + g 







(20) 



resulting in 

d (g) - 2(g) a 
2(d- 1) 



9 = 



Md-l)(g) 2 G 



'1 + 



[d(g) 2(g) a] 



(21) 



The percolating case is easily recovered in the limit 
Ei — » +00, gi — > 0, giving the wellknown result ([14]) 

dp2 - 1 



1 



-92 



P c 2 = 1/d 



(22) 
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while simple calculations can show that in the two- 
dimensional symmetric case, d = 2, p = 1/2, g* equals 
the geometric mean ^Jg~ig~2, as stated in the main text. 
Eq.(f2~Tj) is finally used to calculate the effective diffusion 
constant 



D* = 



D 



Dg* 



d(g) - 2(g) a 
2(d-l)(g) 



(23) 



4(d-l)(g)^ 

[d{g) - 2(g) a] 



Following the strategy outlined in [11], we now introduce 
a free parameter g' and define the quantity 



B(g') 



(d-l)g'+ g/ 
According to [11], the following inequalities hold: 
dB(g') 



(24) 



1 - B(g') 
dB(g' 



if g' < min(g) = g 1 



9* > S'(l 

g* < 4 1 + I^vt ) if - 9 ' >max ^) 

Denoting 

dB( 9l ) 



92 



9*i = 9i{l 
9*2 = .92(1 



1 - B( 9l ] 
dB(g 2 ) 



1 - B(g 2 
we get the following inequalities 



9 > 9\=9\ 

* , * 

9 < 92= 92 
and define 
D\ 

D* 



a 2 9l 
a 2 g* 



D 



giving the bounds 



1 

1 - 

Ml 

{9) 

D— 

(9) 



dgi +p(92 - 9i) 
d P(92 - 9\) 



dg2 - (1 -P){92 ~ 9i) 

d0--p)(g<2 -91) 



dgi +p(92 - gi) 
dp{92 - gi) 



(25) 
(26) 

(27) 
(28) 

(29) 



dg2 - (1 -p)(92 - gi) 



D{ <D* <D* 2 



(30) 



Through a simple rescaling of the units we can fix D = 
<7i = 1 and express <?2 via the contrast x — .92/31 and g* 
via its rescaled value g* — g* /gi, to get a better graphical 
representation. The corresponding boundaries and the 
EMA result for D* then read 



jj* = d[p + (1 - p)x] - {x + 1) 
1 2(d- l)[p + (l-p)x\ 

I + 4/1 



(31) 



A(d- l)x 



[ d \p+(l-p) x }-( x +l)Y 



D\ 

m 



1 



p + (1 — p)x 

x 



d{l-p)(x-l) 



d + p(x- 1) 
dp(x — 1) 



with 



p + (1 — p)x I dx — (1 — p)(x — 1) 



b\<b* < b* 2 



(32) 
(33) 

(34) 



These bounds, together with the effective diffusivity of 
eg. pip, are shown in figure 2 of the main text for two- 
and three-dimensional cases and for p = 1/2. 



Truncated exponential distribution 

Let us now consider the case in which the possible val- 
ues of the site energies E\ are exponentially distributed 
in the interval [E 2l E{\ 



P(Ei) 



-/3Ei 



-0E 2 



-/3E 1 



(35) 



Performing the change of variables it is easy to show that 
the bond conductance is a random variable uniformly 
distributed in the interval [g\ = goe~^ El , g 2 = go e_ ^ B2 ]- 
Eq. ([20|) then gives the following condition to be satisfied 
by the effective conductance 



{d-\)g*+g 2 



-92/dg* _ 



(d-l)g* + gi 



-gi/dg* 



(36) 



or by the effective diffusivity D* 



(d-l)(g)D*+g 2 D 



,-g 2 D/{d{g)D*) 



(d-l)(g)D*+ gi D 



- gi D/(d{g)D') 



(37) 



which can be rewritten through the previous rescaling as 



{d-l)(x + l)D* + 2.x 



-2x/{d(x+l)D") 



= (d-l)(x + l)D* + 2\e- 2/Wx+1)f) "> (38) 



The solution D* of this equation is shown graphically in 
Fig. 3 of the main text together with the two bounds we 
get ready to calculate. 
We consider again 



B(g') = 



9-9 



,{d-l) 9 ' + g 
and the inequalities 

9 * > l-B(q') ) if .9'< min (^)=3i 



(39) 



9* < 9 \ l + Y~B^f]) if 9 ' > max (5ij) = 52 
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In this case we have 



B(g') = 1 



dg' 



In 



92 - gi 



(d-l)g' + g 2 



(d-l)g' + 9l 



(40) 



and setting respectively g' = g\ and g' — g 2 , we can write 



B(gi) 



B(g 2 ) 



dg 



.92 - gi 



■In 



dgi + g 2 - gi 



= l--ln(l + yi ) 

yi v > 



dg 2 



.92 - .91 



■In 



dgi 



dg 2 



Ldg 2 - g 2 + gi 
1 + — In (l - y 2 ) 



2/2 

where the couple of parameters 

.92 - 5i 52 - 5i 



2/i 



c?5i 



2/2 



^52 



(41) 



(42) 



(43) 



has been introduced. At the end the two bounds for the 
conductivity are given by 



5 > 5i = 5i 
5 < 5 2 = .92 
from which we obtain 
n* D 

Dl = V) 91 

n* D 

° 2 = V) 92 



2/i 



ln(l + 2/i) 



1-d 
l + d 



/ y2_ 

Vln(l-lft) 



2/i 



ln(l + |/i) 

52 

ln(l - y 2 ) 



- 1 
+ 1 



and the following bounds 

D* x < D* < D* 2 



(44) 
(45) 

(46) 
(47) 

(48) 



Then, fixing again D = g\ = 1, g 2 = x, g* = g*, we write 



IT = 



D* 



x + 1 
2x 



x + 1 



l-d + 
1-d- 



x - 1 



ln(l + ^i) 

x - 1 



xln(l - ^=i' 



da: 



(49) 
(50) 



with 



D{ < D* < D* 2 



(51) 
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